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Preface 


ruder  Contract  No.  F49G20-C-89-0(W'<.  NTNl-’/NORSAR  is  conducting  re¬ 
search  wit  hi n  a  wide  range  of  su  bjects  relevant  to  seismic  monitoring.  I  lie  einpha 
-is  of  the  researt  It  program  is  on  developing  and  assessing  methods  for  processing 
of  data  recorded  by  networks  of  small-aperture  arrays  and  .'{-com  ponen  I  stations, 
forevents  both  at  regional  and  teleseismic  distances.  In  addition,  more  general 
sei.smological  research  topics  are  addressed. 

Kach  quarterly  technical  report  under  this  contract  presents  one  or  several 
separate  in  vest  igat  ions  addressing  specific  problems  within  tin*  scope  of ' ;  he  s(  ate 
me.ut  of  work.  Summaries  of  the  research  efforts  within  tin-  ptogiam  as  a  whole 
are  given  in  annual  technical  reports. 

"I  his  Scientific  Report.  No.  .'.{  is  the  annual  technical  report  for  the  period  1 
January  -  30  September  1989.  It  contains  five  separate  contributions,  addressing 
the  main  topics  within  the  scope  of  work  for  this  contract. 
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1  SUMMARY 


This  Annual  Technical  Report  describes  the  work  accomplished  under  Contract 
No.  F49620-C-89-0038  during  the  period  1  January  -  30  September  1989.  The 
report  comprises  five  separate  investigations,  addressing  the  main  topics  within 
the  scope  of  work. 

As  part  ol  our  effort  to  enhance  the  performance  of  the  Intelligent  Monitoring 
System,  we  have  embarked  on  a  series  of  projects  to  aid  in  the  understanding 
of  the  behavior  of  the  tectonic  processes  currently  acting  within  Norway  as  well 
as  wave  propagation  effects  for  observed  seismic  phases.  In  paragraph  2.1  ol 
this  report,  new  information  is  conveyed  about  earthquake  faulting  mechanisms 
and  quantification  of  the  structural  information  contained  within  waveform  data 
observed  in  Norway. 

An  analysis  of  the  location  capabilities  of  NOR  ESS  and  ARCESS,  and  of 
the  3-component  stations  within  the  arrays,  leads  to  the  following  conclusions 
concerning  3-component  slowness  solutions  for  P  waves:  (1)  There  is  a.  relatively 
large  scatter  in  the  solutions  for  events  from  the  same  source  region,  and  (2) 
there  are  significant  differences  between  the  solutions  at  the  different  stations.  In 
the  contribution  of  paragraph  2.2  of  this  report,  topographic  effects  have  been 
analyzed  by  comparing  slowness  solutions  for  P  from  explosions  in  three  differ¬ 
ent  source  regions.  Numerical  modelling  demonstrates  that  surface  topography 
can  explain  about  half  of  the  slowness  anomalies  observed  on  three-component 
stations. 

Paragraph  2.3  is  a  status  report  on  work  related  to  local  magnitudes  and 
regional  wave  attenuation.  The  work  comprises  the  establishment  of  a  large 
data  base  of  regional  records,  the  development  of  analysis  tools  and  strategies, 
and  subsequent  data  analysis  aimed  at  derivation  of  regional  wave  attenuation 
relationships  and  the  development  of  a  new  magnitude  scale  for  Norway  and 
adjacent  areas. 

The  paper  of  paragraph  2.4  presents  a  method  for  associating  phase  detections 
from  a  network  of  stations,  which  is  analogous  to  the  conventional  delay-and- 
sum  beamforming  commonly  applied  in  array  detection  processing.  Examples 
of  application  of  this  method  are  given,  based  on  data  recorded  at  the  three 
regional  arrays  (NORESS,  ARCESS,  FINESA)  in  Fennoscandia.  The  problem  of 
continuously  monitoring  the  regional  seismic  noise  field  is  also  addressed,  with  the 
purpose  of  obtaining  a  quantitative  assessment  of  the  upper  limit  of  magnitudes 
of  seismic  events  that  would  go  undetected  by  a  given  network. 

Two  investigations  aimed  at  deriving  legion-specific  knowledge  from  analysis 
of  NORESS  data  are  described  in  paragraph  2.5.  The  first  investigation  is  based 
on  a  study  of  the  complete  NORESS  detection  lists  for  the  period  1985-1988,  and 
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statistics  on  phase  velocities  and  arrival  azimuth  residuals  are  obtained  for  events 
reported  in  the  regional  network  bulletins.  The  second  investigation  is  a  detailed 
study  of  103  events  in  the  western  Norway  area.  The  study  shows  that  arrival 
azimuth  residuals  observed  at  NORESS  for  Pn  and  Sn  phases  vary  systematically 
with  source  region  in  the  western  Norway  area. 


2  SUMMARY  OF  TECHNICAL  FINDINGS  AND 
ACCOMPLISHMENTS 

2.1  Surface  Wave  Modelling  for  Source  and  Structure  Inversions 
Introduction 

As  a  part  of  our  effort  to  enhance  the  performance  of  the  Intelligent  Monitor¬ 
ing  System  we  have  embarked  on  a  series  of  projects  to  aid  in  the  understanding 
of  the  behavior  of  the  tectonic  processes  currently  acting  within  Norway  as  well 
as  wave  propagation  effects  for  observed  seismic  phases.  A  proper  under., land 
mg  of  dynamic  processes  within  the  crust  and  upper  mantle  requires  a  detailed 
knowledge  about  the  occurrence  of  earthquakes  and  tlm  medium  through  which 
the  seismic  waves  propagate. 

Norway  and  its  surrounding  areas,  in  particular  the  northern  continental  mar¬ 
gin,  have  until  recently  been  relatively  poorly  covert'd  in  terms  of  seismic  instru¬ 
mentation.  In  the  last  decade  the  situation  has  been  much  improved,  however, 
through  the  installation  of  new  regional  and  local  networks  and  arrays.  This 
increased  number  cf  Nations  has  resulted  in  more  detailed  delineation  of  seismic¬ 
ity  patterns  as  well  as  improved  focal  mechanism  solutions  and  local  magnitude 
scales.  The  methods  of  matching  synthetic  waveforms  to  observed  data  for  the  re¬ 
trieval  of  source  processes  has  further  increased  our  ability  to  discriminat  e  various 
source  types  from  the  fundamental  observations  of  seismic  data.  This  method- 
olgy  is  strongly  coupled  to  the  determination  of  crustal  structure  for  which  we 
have  also  gained  much  information. 

The  purpose  of  this  contribution  is  to  convey  significant  new  information 
about  earthquake  faulting  mechanisms  and  quantification  of  the  structural  infor¬ 
mation  contained  within  waveform  data  observed  within  Norway  and  along  the 
Norwegian  continental  margin  areas.  Discussion  of  the  results  of  the  modelling 
of  specific  events  is  left  to  a  further  report,  or  can  be  found  in  liungum  ct.  al. 
(1989). 

Source  Mechanism  Determination 

For  an  area  such  as  the  Norwegian  continental  shelf,  with  low-to-intermediate 
seismicity  as  described  above,  the  earlier  (pre  1980)  limitations  in  instrumental 
coverage  also  severely  limited  the  possibilities  for  obtaining  information  about 
crustal  stress  conditions  through  earthquake  focal  mechanisms.  More  and  better 
seismological  stations  as  well  as  improved  methods  for  source  analysis  through 
waveform  modelling  have  now  changed  this  situation  markedly. 

A  major  effort  towards  obtaining  new  focal  mechanism  solutions  has  resulted 
in  28  new  solutions.  Focal  mechanism  solutions  for  several  of  the  earthquakes 
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proved  to  be  quite  difficult  to  obtain  using  the  conventional  first  motion  analysis. 
There  either  were  not  enough  data  to  well  constrain  the  nodal  planes  or  errors 
were  introduced  in  the  take-off  angles  from  the  source  due  to  misinterpretation  of 
the  travel  path  of  the  seismic  phase  because  of  the  complicated  geometry  of  the 
crust-mantle  interface  (as  was  the  case  for  the  event  in  Figure  2.1.1,  Event  5  in 
Table  2.1.1).  In  six  of  these  cases,  source  mechanisms  have  been  retrieved  through 
waveform  modelling,  using  the  same  method  as  Hansen  et.  al.  (1989)  applied  for 
two  of  the  largest  earthquakes  in  the  area  (Events  2  and  5  in  Table  2.1.1).  The 
results  have  been  found  to  be  very  stable  even  when  data  from  only  one  station 
have  been  used  (NORSAR  long-period  for  Event  1,  NORESS  intermediate-period 
for  Events  3  and  4). 

Since  several  events  were  found  to  have  produced  fairly  simple  and  well 
recorded  long  period  surface  waves,  an  approach  of  source  mechanism  retrieval 
was  adopted  that  combines  broad-band  waveform  modelling  with  the  more  con¬ 
ventional  first  motion  analysis.  The  method  (Hansen,  1989)  consists  of  low  pass 
filtering  the  broad-band  records  from  the  NORESS  array  (NRS)  and  Kongs- 
berg  (KONO)  to  emphasize  the  low  frequency  waves  up  to  about  6  seconds  pe¬ 
riod  and  then  matching  the  waveforms  to  synthetic  seismograms  computed  for  a 
given  earth  structure  and  source  mechanism  in  both  a  forward  and  inverse  sense. 
Green’s  functions  vere  computed  using  the  Locked  Mode  Approximation  method 
of  Harvey  (1981)  for  a  crust-mantle  structure  derived  from  Mykkeltveit  (1980) 
for  the  crust  and  Stuart  (1978)  for  the  upper  mantle. 

The  results  of  this  modelling  are  illustrated  in  Figures  2.1.1  (forward  mod¬ 
elling)  and  2.1.2  (inverse  modelling)  for  Events  5  and  2.  respectively.  Three- 
component  seismograms  recorded  for  each  event  are  shown  after  filtering  and 
rotation  to  vertical,  radial,  and  transverse  components.  The  corresponding  syn¬ 
thetic  seismograms  are  shown  just  below  the  observed  ones.  The  method  was 
first  applied  to  Event  5  (Figure  2.1.1)  where  excellent  results  were  obtained.  The 
inverse  procedures  were  then  developed  and  the  data  set  for  the  event  in  Fig¬ 
ure  2.1.2  was  augmented  to  include  the  three  component  Kongsberg  (KONO) 
long  period  data.  To  the  right  of  Figures  2.1.1  and  2.1.2  are  the  focal  mecha¬ 
nism  drawings  (lower  hemisphere  stereographic  projections)  for  the  nodal  planes 
obtained  from  the  waveform  modelling  together  with  the  first  motion  data.  It 
is  obvious  that  the  first  motion  data  are  to  some  extent  discordant,  either  due 
to  incorrect  readings  or  false  interpretations.  Experimenting  with  different  ve¬ 
locity  models  for  location,  and  varying  the  source  depth  failed  to  improve  the 
interpretation  of  the  first  motion  data.  However,  they  were  useful  in  helping  to 
discriminate  between  say  two  different  waveform  solutions  and  serve  to  supply 
constraints  on  further  interpretation  of  velocity  models  in  the  area.  Several  of 
these  events  were  verified  with  forward  waveform  modelling  even  if  the  inversion 
procedure  was  not  deemed  warranted  by  the  data.  As  an  example,  Figures  2.1.3 
and  2.1.4  illustrate  the  fit  from  forward  modelling  of  two  earthquakes  (Events  4 
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and  1  in  Table  2.1.1)  near  the  Lofoten  Islands  at  a  distance  of  more  than  800 
kilometers  from  the  recording  stations  at  NORESS  (Figure  2.1.3)  and  NORSAR 
(Figure  2.1.4).  It  is  very  encouraging  to  find  that  initial  interpretations  of  the 
sparse  amount  of  first  motion  data  can  be  verified  to  a  much  higher  degree  of 
confidence. 

Due  to  the  complexity  of  the  crust  and  upper  mantle  structure  in  the  vicinity 
of  the  earthquakes,  it  was  necessary  to  explore  the  effect  of  the  assumed  veloc¬ 
ity  structure  on  the  low  frequency  waveform  modelling.  Since  we  are  observing 
phases  with  wavelengths  on  the  order  of  20  kilometers  and  longer,  much  of  the 
complexity  affecting  the  high  frequency  first  motion  data  is  smoothed  out.  The 
structures  used  for  computing  synth-  tie  sol  mograms  therefore  reflected  the  aver 
age,  and  more  homogeneous,  prop*  rties  sampled  along  the  majority  of  the  travel 
path  from  the  earthquake  to  the  receiving  stations.  By  observing  the  change  in 
the  synthetic  waveforms  as  a  function  of  structural  model  and  earthquake  sottne 
depth  it  was  verified  that  the  solutions  for  the  focal  mechanisms  and  source  depths 
are  quite  robust  for  the  ’ong  period  data.  Variations  in  the  compression  and  ten¬ 
sion  axes  of  the  focal  sphere  were  found  to  be  less  than  about  5  degrees  du«*  to 
change  .  in  velocity  models  for  a  single  station  solution,  and  wore  improved  with 
the  inclusion  of  a  second  three-component  station  (KONO).  It  was  also  found 
that  source  depths  could  be  verified  from  the  modelling  to  a  precision  of  about 
±  5  km. 

For  Event  2  (Fig-  2.1.2),  the  waveform  fit  is  excellent.  It  is  seen  from  the  focal 
mechanism  plot  that  the  fault  plane  solution  would  be  difficult  to  obtain  from  first 
motion  data  alo  t.  In  fact,  an  earlier  interpretation  of  the  first  motion  data  of  this 
event  indicated  a  normal  fault  with  a  near  vertical  axis  of  maximum  compression. 
However,  the  waveform  modelling  shown  here  for  this  event  completely  rules  out 
this  typo  of  solution.  Instead,  the  solution  for  this  earthquake  dearly  indicates 
oblique  thrust  faulting  in  response  to  a  regional  stress  field  exhibiting  NW-SE 
horizontal  compression. 

For  most  of  the  28  new  events,  a  similar  amplitude  modelling  has  not  been 
feasible  because  of  the  very  different  character  of  the  long  period  waveforms. 
Small  events  did  not  have  enough  energy  in  the  long  periods  for  them  to  be  ob¬ 
served  above  the  noise  levels,  while  some  larger  events  produced  far  too  complex 
waveforms  to  be  easily  interpreted.  The  complex  seismograms  are  likely  a  result 
of  the  rapidly  changing  crustal  structure  between  events  out  to  sea  and  th  -  seis¬ 
mometers  located  in  southeastern  Norway.  In  these  cases,  a  solution  using  first 
motions  only  is  given,  where  a  combination  of  local  and  teleseisrnic  data  helps 
in  constraining  the  nodal  planes.  However,  when  the  events  are  within  local  dis¬ 
tance  ranges,  the  higher  order  surface  wave  trains  can  be  summed  to  produce 
body- waves  that  can  be  fit  to  the  observations  up  to  2  or  3  Hertz.  Figure  2.1.5 
shows  an  example  of  such  an  event  (Event  6)  occurring  very  near  the  NORESS 
and  NORSAR  arrays  with  recording  stations  ranging  from  about  15  to  85  kilome- 
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ters.  Figure  2.1.5  shows  that  the  velocity  model  is  only  crucial  to  events  at  large 
distances  at  higher  frequencies.  An  approach  to  improve  this  result  is  outlined  in 
the  next  section  as  an  improvement  to  the  velocity  models. 

Improved  Velocity  Models 

Although  we  have  achieved  good  results  in  inverting  for  source  parameters  by 
matching  synthetic  an  1  real  seismograms,  our  success  has  been  limited  since  we 
have  only  been  able  to  do  this  for  a  small  number  of  events  and  at  low  freq  encies. 
Since  we  are  using  a  laterally  homogeneous  modelling  technique,  it  can  be  argued 
that  such  a  technique  may  never  be  suitable  for  inversion  in  situations  where 
the  structure  depth  dependence  changes  significantly  aiong  the  source-receiver 
propagation  path,  such  as  earthquakes  in  the  oceanic  crust  that  are  recorded 
within  the  continental  crust.  Unfortunately  the  p  oblem  of  regional  forward 
wave  propagation  modelling  through  arbitrary  inhomogeneous  media  is  effectively 
intractable  at  intermediate  and  high  frequencies,  and  even  at  low  frequencies  it 
can  only  be  accomplished  by  making  assumptions  on  the  nature  of  the  media 
inhomogeneities,  such  as  slow  and  weak  lateral  variations,  that  are  often  violated 
in  real  world  situations. 

The  success  that  we  have  had  demonstrates  that,  at  least  in  certain  situations, 
laterally  homogeneous  modelling  can  be  used  to  “image  out”  both  the  structure 
and  the  source  parameters.  These  situations  seem  to  involve  cases  where  ti  e 
source- recei  ver  propagation  paths  stay  within  the  boundaries  of  the  continental 
crust,  which  lends  more  credence  to  the  laterally  homogeneous  assumption.  There 
are  not  many  crustal  events,  however,  that  are  large  enough  to  produce  significant 
energy  in  the  low  frequencies.  We  would  like  to  be  able  to  investigate  small 
crustal  events,  such  as  local  and  near-regional  explosions,  but  we  can  only  do  thi 
by  looking  at  higher  frequencies  where  the  source  produces  significant  energy. 
When  we  look  at  th  •  events  that  we  can  match  well  at  low  frequencies,  our 
matches  deteriorate  rapidly  as  we  increase  the  frequency  bandwidth,  which  is 
not  surprising  considering  the  relatively  crude  nature  of  the  structural  models 
that  we  have  been  using.  As  we  raise  the  frequency  content  of  our  synthetic 
seismograms,  we  find  that  they  become  increasingly  sensitive  to  detailed  changes 
in  the  structural  model.  Thus  it  becomes  <  .ident  that,  in  order  *o  match  synthetic 
and  real  data  at  higher  frequencies,  it  will  be  necessary  to  invert  for  structural 
parameters  as  well  as  for  source  parameters. 

We  have  begun  a  systematic  approach  to  the  comparison-inversion  problem 
which  will  take  advantage  of  some  recent  work  that  has  been  done  in  this  area 
and  winch  should  provide  a  clear  and  realizable  path  to  detailed  s'ructural  mod¬ 
els  that  will  produce  optimal  fits  between  the  real  and  synthetic  seismograms. 
The  core  of  our  me* hod  involves  a  modified  version  of  the  direct  inversion  tech¬ 
nique  developed  by  Gomberg  and  Masters  (1988)  which  uses  iterative  linear  in¬ 
version  based  upon  differential  seismograms  that  are  analytically  computed  from 
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a  Locked  Mode  representation  of  the  solution.  We  will  carry  out  this  inversion 
in  stages,  starting  with  very  low  frequency,  band  limited  comparisons  of  real  and 
synthetic  seismograms  which  will  determine  the  gross  characteristics  of  the  struc¬ 
ture.  We  will  then  proceed  to  open  up  the  bandwidth  of  comparison  in  steps  and 
invert  for  increasingly  detailed  aspects  of  the  structural  model  while  maintaining 
the  average  characteristics  that  were  determined  in  the  lower  band  inversions  for 
both  source  and  structure.  At  the  same  time  it  may  be  necessary  as  we  go  u p 
in  frequency  to  re-estimate  the  source  parameters  for  each  more  complex  struc¬ 
ture  iteration.  In  this  manner  we  hope  to  obtain  the  best  average  velocity-depth 
dependencies  in  specific  tectonic  provinces. 

V\e  have  implemented  the  core  computations  for  the  differential  seismograms 
at  XORSAR  and  (.'IRES  based  on  the  Locked  Modi'  methods  of  llarvev  ( 1  OS  1 ) . 
At  this  stage'  we  are  abb'  to  compute  perturbation  seismograms  that  can  be  used 
to  determine  the  effect,  on  the  seismogram  of  small  changes  in  layer  velocities. 
I'he  current  implementation  is  functionally  identical  to  that  of  Gomberg  and 
Masters  in  that  it  is  a  perturbation  implementation  which  assumes  that  only 
the  eigenvalues  change,  as  opposed  to  a  true  differential  implementation  which 
accounts  for  both  changes  in  the  eigenfunctions  and  eigenvalues.  The  pertur¬ 
bation  seismograms  become  inaccurate  as  the  source- receiver  distance  decreases 
and  the  frequency  increases,  so  we  are  in  the  process  of  implementing  the  exact 
differential  seismogram  computations  and  undertaking  a  detailed  study  to  show 
where  the  perturbation  approximation  breaks  down.  We  have  implemented  the 
core  differential  computations  in  a  significantly  different  manner  from  that  of 
Gomberg  and  Masters  which  allows  us  to  do  the  exact  differential  computations 
in  a  straightforward  manner. 


Roger  A.  Hansen 

Danny  J.  Harvey,  (’IRES,  University  of  Colorado 
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No. 

Reg. 

Date 

Lat 

N 

Lon 

E 

Depth 

Mag. 

T- 

Az 

axis 

Dip 

P- 

Az 

axis 

Dip 

Type 

1 

NN 

810903 

69.62 

13.68 

1 2C 

4.7 

49 

38 

3IM 

22 

RO 

2 

WN 

86020") 

62.71 

4.69 

20C 

4.9 

198 

55 

303 

10 

R 

3 

VVN 

861026 

61.83 

3.20 

14B 

4.5 

305 

78 

112 

12 

H 

1 

NN 

880131 

6X. 03 

9.58 

20C 

4.3 

92 

24 

329 

51 

NO 

5 

WN 

890123 

61.97 

4.42 

26  A 

5  1 

239 

72 

124 

8 

R 

6 

SN 

890110 

60.61 

11.40 

22B 

1.9 

35 

42 

303 

3 

RO 

Table  2.1.1.  Earthquake  focal  mechanisms  collected  for  this  study,  with  the 
different  columns  indicating  event  number,  region  (SN,  Southeastern  Norway. 
VY.N,  Western  Norway,  NN,  Northern  Norway),  date  (year,  month,  day),  latitude, 
longitude,  focal  depth  (A,  precision  <  3  km;  B,  4-6  km;  C,  7-9  km,  1),  >  10 
km),  magnitude  (M^),  T-axis  (minimum  compression)  azimuth  and  dip,  E-axis 
(maximum  compression)  azimuth  and  dip,  type  of  mechanism  (  N,  normal;  NO, 
normal  oblique;  R,  reverse  or  thrust;  RO,  reverse  oblique). 
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Date:  01  23  89  1406 


Time,  (sec) 


Filtering  0  0.00  8  0.07 

Green  Function  source  depth  =  30.0 

nlvzlhz30jan23.nofilt.grn  9/22/8918:58:37  N»8»A5 


Figure  2,1.1.  Focal  mechanism  solution  from  forward  waveform  modelling, 
F,vont  5,  2,'{  January  1989.  The  left  frame  shows  three  components  (vortical,  ra 
dial  and  transverse)  of  NOKISSS  broad  hand  data  together  with  theoretical  seis¬ 
mograms  from  the  double  couple  moment  tensor  immediately  below  each  trace, 
while  the  right  frame  shows  the  double  couple  solution  from  the  moment  tensor. 
T  he  first  motion  data  not  read  by  the  authors  are  indicated  by  plus  and  minus 
signs  for  compressions  and  dilations,  respectively. 


10 


Date:  02  05  86  1753 


Time,  (sec) 


Filtering  0  0.00  2  0.07 
norsar  1 hz25feb5kono.grn 


3/22/8918:17:44  N08SAS 


Figure  2.1.2.  Focal  mechanism  solution  from  inversion  of  waveform  data, 
Event  2,  5  February  1986.  The  left  frame  shows  three  components  (vertical, 
radial  and  transverse)  of  NORESS  broad  band  data  and  Kongsberg  Long  Period 
data  together  with  theoretical  seismograms  from  the  deviatoric  moment  tensor 
immediately  below  each  trace,  while  the  right  frame  shows  the  best  double  couple 
solution  from  the  moment  tensor.  The  percent  non-double  couple  for  this  event 
at  these  low  frequencies  is  less  than  1. 
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Time,  (sec) 


Filtering  0  0.00  4  0  07 

Creen  Function  source  depth  =  20.0 

nlvzlhz20jan31  grn  3/22/8917:34:24  NW8&A3 


Figure  2.1.3.  Focal  mechanism  solution  from  forward  waveform  modelling  for 
a  Lofoten  event  of  31  January  1988.  This  figure  illustrates  the  application  of 
the  modelling  technique  as  applied  to  the  NORESS  Broad-Band  recording  of  an 
event  in  Northern  Norway  at  a  distance  of  more  than  800  km  in  the  magnitude 
range  of  4.0  to  4.5.  The  left  frame  shows  three  components  (vertical,  radial 
and  transverse)  of  NORESS  filtered  Broad-Band  data  together  with  theoretical 
seismograms  immediately  below  each  trace,  while  the  right  frame  shows  the  focal 
mechanism  solution  used  in  the  waveform  modelling. 
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Date:  09  03  81  1839 


Time ,  (sec) 


Fillermg  0  0.00  4  0.10 

Creen  Function  source  depth  =  10  0 

nivzl  hzl0scp3.grn  10/17/8911:34:40 


Figure  2.1.4.  Focal  mechanism  solution  from  forward  waveform  modelling  for 
a  Lofoten  event  of  3  September  1981.  This  figure  illustrates  the  application  of 
the  modelling  technique  as  applied  to  the  NORSAR  LP  recording  of  an  event  in 
Northern  Norway  at  a  distance  of  more  than  800  km  in  the  magnitude  range  of 
4.5  to  5.  The  left  frame  shows  three  components  (vertical,  radial  and  transverse) 
of  NORSAR  Long  Period  data  together  with  theoretical  seismograms  immedi¬ 
ately  below  each  trace,  while  the  right  frame  shows  the  focal  mechanism  solution 
used  in  the  waveform  modelling.  This  figure  illustrates  the  dominance  of  the  Love 
wave  component  for  solutions  derived  from  waveform  modelling. 
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Time  (sec) 


Filtering  0  0  00  6  2.00 

Green  Function  source  depth  =  22.0 

nnla.22  gund  grn  1  I  / 1  0/89 1 4: 1  0:47  W<95g.A8 


Figure  2.1.5.  Focal  mechanism  solution  from  forward  waveform  modelling  for 
a  magnitude  2.1  event  located  within  the  NORSAR  array.  This  figure  illustrates 
the  application  of  the  modelling  technique  as  applied  to  the  NORSAR  recording 
of  an  event  at  a  distance  of  40  km  and  to  a  frequency  of  3  hz.  The  left  frame 
shows  the  vertical  component  of  NORSAR  short  period  data  together  with  a 
synthetic  seismogram  calculated  with  the  fault  plane  solution  from  the  right  frame 
constrained  by  the  first  motion  data. 


11 


2.2  Topographic  effects  on  arrays  and  three-component  stations 

An  analysis  of  the  location  capabilities  of  NORESS  and  AIK'ESS.  and  of  the 
3-componcnt  stations  within  the  arrays,  leads  to  the  following  conclusions  con¬ 
cerning  3-component  slowness  solutions  for  P  waves:  (1)  There  is  a  relatively 
large  scatter  in  the  solutions  for  events  from  the  same  source  region,  and  (2) 
there  are  significant  differences  between  the  solutions  at  the  different  stations.  A 
third  result  is  common  for  both  array-  and  3-component  solutions:  (3)  There  are 
signicant  anomalies  in  the  slowness  of  P  from  particular  source  regions.  For  a 
proper  evaluation  of  NORESS  and  similar  arrays,  it  is  important  to  understand 
the  cause  of  these  observational  phenomena.  In  this  work  we  have  analyzed 
topographic  effects  by  comparing  slowness  solutions  for  P  from  explosions  in  3 
different  source  regions:  near  Leningrad,  in  E.  Kazakh,  and  near  Rlasjo  in  S.W. 
N  orway. 

Slowness  analysis 

A  unified  approach  to  slowness  analysis  with  arrays  and  3-component  stations 
is  possible  by  expressing  solutions  in  terms  of  a  covariance  matrix  C  (e.g.,  Es- 
mersoy  et  al ,  1985).  Here  we  introduce  C  as  a  function  of  slowness  »  by  phase 
shifting  the  signals: 

Cnm(s)  =  J  Fn(xv,s)F^(w,s)dw/2TT  (1) 

where 

Fn(w,s)  =  Fn{w)  exp(iws.xn) 

and  Fn  is  the  Fourier  spectrum  at  channel  n.  Using  C  of  eq.  (1),  the  generaliza¬ 
tion  of  conventional  beamforming  is  given  by  the  normalized  response 

P(s)  =  g^g/d  g  |2  frC}  (2) 

where  g  is  the  predicted  displacement  vector  for  slowness  s.  The  generalization 
(2)  can  be  interpreted  as  a  matched  filter  since  the  response  P  depends  on  our 
choice  of  matching  the  covariance  matrix  C: 

gfCg  =  tr ( C  G)  with  G  =  gg1 

Thus  for  a  1-component  array:  gT  =  (1,...,1)T,  and  for  a  3-component 
sensor:  gr  =  (gz,gy,gz)  =  displacement  vector.  The  latter  is  a  function  of 
slowness  s,  and  the  surface  interaction  must  be  taken  into  account.  On  this 
account  we  may  expect  3-component  slowness  solutions  to  be  relatively  sensitive 
to  the  choice  of  near-surface  model.  It  is  possible  to  extend  the  procedure  to  a 
generalization  of  optimum  beamforming.  For  example,  the  normalized  response 
of  the  maximum  likelihood  method  is 

P'(s)  =  {^C-'g}-’  |  g  |2  /<rC  (3) 
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Although  eq.  (3)  leads  to  solutions  with  apparently  higher  resolution,  the 
solutions  are  less  stable  and  for  location  purposes  stability  is  more  important. 
Hence  we  will  proceed  the  analysis  based  on  eq.  (2). 

The  Leningrad  events 

Figure  2.2.1  shows  a  typical  3-component  record  at  NORESS  from  a  mining 
explosion  near  Leningrad.  Five  such  events  were  recorded  with  good  signal-to 
noise  ratio,  and  Figure  2.2.2  summarizes  both  the  array  slowness  solutions  and 
the  3-coinponent  solutions  for  4  stations  within  NORESS.  Standard  error  bars 
indicate  the  scatter  of  solutions  for  the  5  events.  Error  bars  of  3-component 
solutions  are  given  for  only  one  station,  but  the  scatter  for  the  other  stations  is 
similar.  The  theoretical  results  included  in  Figure  2.2.2  will  be  discussed  later. 
The  3-component  slownesses  show  a  relatively  large  scatter,  both  of  solutions  at 
different  stations  for  the  same  event  and  of  solutions  for  different  events  at  the 
same  station.  These  solutions  take  into  account  wave  interaction  with  a  plane 
surface.  An  obvious  extension  is  to  include  the  effect  of  surface  topography. 

Interaction  with  surface  topography 

The  usual  correction  for  surface  topography  implies  a  time  correction  for  ele¬ 
vation,  and  possibly  a  particle  motion  correction  for  surface  slope.  These  correc¬ 
tions  are  consistent  with  geometrical  ray  theory,  which  requires  that  topographic 
relief  be  smooth  on  the  scale  of  a  wavelength.  However,  the  topography  of  the 
NORESS  area  is  not  smooth  on  the  scale  of  the  relevant  wavelengths  (~  2  km), 
and  wave  scattering  effects  may  be  important.  The  same  is  true  for  the  ARCESS 
area.  To  evaluate  the  scattering,  we  have  applied  a  perturbation  method  where 
the  solution  is  obtained  recursively  in  wavenumber  space  (Doornbos,  1988).  For 
the  present  purpose  we  need  the  surface  displacement  u(x,y).  The  required  per¬ 
turbation  series  is  then 

GO 

U(fc,,*y)=E  U<n>(fcr,fc„)  (4) 

n=0 

where  U(A:X,A.'V)  is  the  Fourier  transform  of  u(x,y).  and  U,ri^  is  a  function  of 
U"~m,  1  <  m  <  n,  of  the  surface  topography /(i,  y),  and  of  the  surface  gradients 
df/dx,  df/dy.  The  zeroth  order  term  U°  gives  the  conventional  free  surface 
response  for  a  plane,  t)ie  first  order  term  U1  includes  the  Born  approximation, 
and  the  higher  order  terms  account  for  multiple  scattering. 

We  have  synthesized  the  free  surface  response  at  NORESS  in  the  signal  band¬ 
width  (2-4  Hz),  assuming  the  incident  wave  direction  is  given  by  the  array  slow¬ 
ness  solution.  We  then  applied  the  method  of  equation  (2)  to  the  3  components  of 
surface  displacement,  to  obtain  theoretical  slowness  solutions  for  the  3-component 
station  sites.  The  results  are  included  in  Figure  2.2.2.  These  results  explain  about 
half  of  the  observed  anomalies.  We  speculate  that  shallow  subsurface  structure 
may  enhance  the  surface  topographic  effect.  The  results  represent  a  weighted 
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integral  over  the  bandwidth  2-4  Hz.  but  single  frequency  solutions  give  a  signifi¬ 
cant  variation  with  frequency.  This  is  consistent  with  the  notion  of  response  by  a 
rough  surface  as  an  interference  pattern.  Thus,  different  signal  spectra  from  the 
same  source  region  may  lead  to  different  slowness  solutions,  which  might  explain 
the  scatter  of  solutions  at  a  single  station. 

E.  Kazakh  events 

Teleseismic  P  waves  from  recent  nuclear  explosions  in  E.  Kazakh  have  been 
recorded  both  at  NOR  ESS  and  ARCESS.  The  P  wave  spectra  from  all  events  are 
similar  at  NOR  ESS,  but  at  ARCESS  we  can  distinguish  two  groups  of  events.  The 
spectra  within  a  group  are  similar,  but  there  are  significant  differences  between 
the  two  groups.  Representative  spectra  at  ARCESS  are  shown  in  Figure  2.2.;!, 
and  average  slowness  solutions  for  each  group  arc  plotted  in  Figure  2.2.4  (a  and 
b).  There  is  a  slight  difference  between  the  array  slowness  solutions,  but  there  is 
a  large  difference  between  the  3-component  solutions  for  the  two  groups.  \Ve  have 
also  plotted  the  theoretical  slowness  solutions  for  the  3-component  station  sites. 
The  assumed  signal  bandwidth  was  0.9-2. 5  Hz  for  the  ‘‘high-frequency”  group  of 
events,  and  0.9- 1.6  Hz  for  the  “low-frequency”  group.  The  results  suggest  that 
the  surface  topography  at  ARCESS  explains  slightly  less  than  half  of  the  observed 
anomalies.  It  is  also  clear  that  the  slowness  results  for  the  two  spectra  are  rather 
different  in  accordance  with  observations.  This  supports  the  explanation  that 
the  often  significant  variation  of  3-component  slowness  solutions  for  signals  from 
events  in  the  same  source  region  can  be  related  to  differences  in  the  signal  spectra. 

Blasjo  events 

We  have  analyzed  in  some  detail  the  NORESS  records  of  Pn  from  a  suite  of 
6  mining  explosions  in  the  Blasjo  area  in  S.W.  Norway.  The  mining  site  is  about 
300  km  from  NORESS  in  an  azimuth  direction  of  about  240°.  A  record  section 
for  one  of  these  events  with  a  plot  of  slowness  as  a  function  of  time  is  shown  in 
Figure  2.2.5.  Slowness  and  azimuth  of  the  first  arrival  are  consistent  with  Pn  in  a 
one-dimensional  crust-mantle  model  (e.g.,  Menke  and  Richards,  1980).  However, 
it  can  be  seen  that  the  first  arrival  is  relatively  weak.  In  fact,  this  arrival  is  easily 
missed  for  small  events,  whence  Pn  is  often  associated  with  the  dominant  part  of 
the  wave  train.  The  dominant  part  about  0. 5-0.6  seconds  after  the  first  arrival 
has  a  consistent  slowness  and  azimuth  anomaly  for  all  events.  Slowness  solutions 
as  a  function  of  frequency  as  summarized  in  the  slowness/azimuth  spectrum  of 
Figure  2.2.6  show  that  the  anomaly  is  related  to  the  frequency  range  2-4  Hz; 
this  is  also  the  range  where  the  signal  has  its  maximum  energy.  The  nearly 
plane  wavefront  indicates  that  the  anomaly  cannot  be  generated  near  the  surface, 
and  it  is  suggested  that  these  “Pn”  waves  are  actually  the  result  of  scattering 
at  depth.  From  the  measured  slowness  and  time  delay,  the  scattering  source 
is  constrained  to  be  within  the  depth  range  of  the  Moho  (30-40  km).  Hence 
the  observational  results  are  consistent  with  scattering  by  topographic  relief  of 
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the  Moho.  An  interesting  geological  aspect  is  that  the  inferred  location  of  the 
proposed  topographic  feature  coincides  with  the  border  of  the  Oslo  Graben. 

Scattering  by  Moho  topography 

Scattering  will  of  course  afToct  all  waves  interacting  with  a  rough  Moho  dis¬ 
continuity.  However,  the  scattered  waves  usually  arrive  in  the  coda  of  a  relatively 
strong  primary  wave.  In  contrast,  the  first  arriving  Pn  is  relatively  weak  du<' 
to  the  small  coefficient  of  refraction  through  the  Moho,  and  scattering  due  to 
topography  of  the  boundary  may  dominate  the  wave  train. 

To  illustrate  these  concepts,  we  have  calculated  generalized  transmission  coef¬ 
ficients  for  a  rough  Moho.  In  analogy  to  equation  (4)  for  the  displacement  vector, 
we  introduce  a  recursive  solution  for  the  vector  of  scattering  coefficients: 

B(kr,ky)=  jr  (5) 

Ti  —  0 


The  components  of  B*,ud  are  just  the  plane  wave  transmission  coefficients 
for  a  plane  interface.  From  the  components  of  B  we  can  calculate  the  energy 
flux  for  any  wave  type.  Figure  2.2.7  shows  the  energy  flux  of  P  transmitted 
upward  through  a  rough  Moho,  as  a  function  of  incident  P  below  the  boundary. 
The  topography  here  is  characterized  by  a  correlation  length  of  5.6  km  and  an 
average  height  of  1  km,  and  the  wave  frequency  is  3  Hz.  Two  modes  of  scattering 
are  shown:  (1)  The  specular  flux  E°  in  the  direction  defined  by  the  plane  wave- 
plane  interface  concept.  The  specular  flux  through  a  rough  interface  is  reduced 
with  respect  to  the  flux  through  a  plane.  (2)  The  diffuse  flux  E,c  due  to  multiple 
scattering  in  all  upward  directions.  E,c  does  not  exist  for  a  plane  interface. 
The  figure  illustrates  well  the  sharp  increase  of  the  ratio  Esc / E°  as  the  slowness 
approaches  the  critical  value  corresponding  to  Pn,  thus  supporting  in  a  qualitative 
way  the  scattering  model  for  propagation  of  this  wave.  The  results  suggest  that 
a  careful  calibration  is  needed  before  using  this  phase  for  event  location  and 
velocity  determination  purposes. 

D. J.  Doornbos 

E.  0degaard,  Univ.  of  Oslo 
T.  Kvaerna 
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Fig.  2.2.1.  1  ypical  three-roiiiponont  NORESS  records  from  event  in  Leningrad 
region.  Scaling  factors  of  different  components  are  shown  to  the  '  >ft. 
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Fig.  2.2.2.  Slowness  solutions  at  NORESS.  Average  over  P  from  5  events  near 
Leningrad.  Array  solution  (ALLV)  with  standard  error  bars.  Three-romponent 
solutions  at  the  indicated  sites,  with  standard  error  bars  for  site  C2.  Theoretical 
solutions  including  response  to  surface  topography  are  framed. 


21 


Sy  (s/km)  Sy  (s/km) 


0.02 


a) 


•  -  ALLV 
V  -  A03C 

♦  -  C23C 
a  -  C41C 

*  -  C'  7  3  C 


0.01 


0.00 


0.07  0 . 08  0.09  0^0  0.11 


Sx  (s/km) 


0.02 


#  -  ALLV 
V  -  A03C 

♦  -  C23C 
-  C43C 

a  -  C73C 


0.01 


0 


0 . 00 


-0.01 


I 

i 


-0.02 

0 . 07 


0 . 08 


j  —  | — 

0-09  0.10 

Sx  (s/km) 


Fig.  2.2.4.  Slowness  solutions  at  ARCESS.  Details  as  in  Fig.  2.2.2. 

a.  Average  over  4  events  from  group  A  in  Kazakh. 

b.  Average  over  5  events  from  group  B  in  Kazakh. 
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Fig.  2.2.5.  NORESS  record  of  Pn  from  an  explosion  in  the  Blasjo  area.  The  top 
trace  is  a  single  channel  seismogram.  Array  solutions  of  slowness  and  azimuth 
as  a  function  of  time  are  shown  below.  The  symbol  size  in  the  plot  indicates  the 
relative  signal  power. 
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Fig.  2.2.6.  Slowness  and  azimuth  solutions  of  Pn  as  a  function  of  frequency, 
from  6  explosions  in  the  Blasj®  area.  The  symbol  size  indicates  the  relative  signal 
power. 
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Fig.  2.2.7.  Energy  flux  at  3  Hz  through  rough  solid-solid  interface  with  the 
velocity-density  structure  v^~  =  6.8/8. 1  km/s,  -  3. 7/4.5  km/s,  p+/~  = 
2.9/3. 4  g/cm3.  The  boundary  roughness  is  characterized  by  an  average  height  of 
1  km  and  a  correlation  length  of  5.6  km.  The  scattered  flux  is  upward  through 
the  boundary. 

-  :  specular  direction; - :  integrated  diffuse  flux. 


2.3  Local  magnitudes  and  regional  wave  attenuation 

Introduction 

The  research  conducted  during  the  last  year  within  the  field  of  regional  wave 
attenuation  and  local  magnitudes  has  been  organized  as  part  of  a  more  long  term 
plan,  with  the  following  done  so  far: 


•  Establishment  of  a  large  data  base  of  regional  records,  stored  in  a  common 
format 

•  Development  of  analysis  tools  and  strategics,  and  preliminary  analysis  of 
data 


In  the  following,  we  will  describe  each  of  these  subjects  in  more  details. 

Data  base  of  records 

For  a  long  time,  the  network  of  seismic  stations  in  Norway  was  limited  to 
a  few  (2-4)  conventional  analog  stations,  and  from  1970,  the  NORSAR  array. 
Then,  between  1980  and  1983,  the  regional  Southern  Norway  Seismic  Network 
(SNSN)  was  in  operation,  the  Western  Norway  Seismic  Network  (WNSN)  started 
recording  in  late  1984,  and  the  Northern  Norway  Seismic  Network  (SEISNOR) 
in  early  1987.  With  NORESS  (1985)  and  ARCESS  (1987)  in  addition,  the  whole 
country  and  surrounding  areas  are  therefore  now  quite  well  covered. 

It  is  being  recognized  now,  more  and  more,  that  well-calibrated  magnitude 
scales  are  quite  important  in  very  many  seismological  problems  and  applications. 
This  is  obviously  true  for  problems  related  to  detectability  and  yield  estimation, 
but  also  to  wave  attenuation  where  magnitude  in  some  cases  is  an  important 
scaling  parameter.  It  is  also  being  recognized  that  the  development  of  a  new 
regional  magnitude  scale  requires  a  very  thorough  development  of  its  wave  at¬ 
tenuation  terms,  and  that  this  requires  large  amounts  of  observed  data  in  order 
to  be  sufficiently  reliable.  The  data  situation  described  above  for  Norway  makes 
this  the  right  time  for  pursuing  these  goals. 

To  this  end,  we  have  undertaken  a  major  effort  in  terms  of  establishing  a  large 
data  base  of  regional  earthquake  records,  later  to  be  used  in  the  investigation  of 
a  variety  of  seismological  problems.  The  data  collected  so  far  consists  of  more 
than  200  earthquakes  and  more  than  1000  individual  records,  with  magnitudes 
in  the  range  Mi  1. 7-5.3  and  epicentral  distances  in  the  range  10-2000  km.  All  of 
these  records  have  been  converted  to  a  common  format,  convenient  for  interactive 
analysis  by  program  packages  which  already  include  all  necessary  system  response 
information. 
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For  the  records  selected,  all  epicenters,  recording  stations  and  ray  paths  are 
shown  in  Fig.  2.3.1,  while  the  magnitude-distance  distribution  is  shown  in  Fig. 
2.3.2.  It  should  be  noted  here  that  some  of  these  records  (near-field,  large  mag 
nitudes)  are  clipped,  while  others  (far  field,  low  magnitudes)  may  be  too  much 
affected  by  noise. 

Analysis  strategies 

In  following  Richter  (1935;  1958),  the  local  magnitude  Mi  is  defined  as  (see 
also  Boore,  1989): 

Mi  =  log  A  -  logAo  +  S  (1 ) 

where  A  is  measured  amplitude  (0-p)  in  mm  on  a  Wood-Anderson  seismometer 
recording,  S  is  a  station  correction,  and 

-  logAo  =  a  ■  fop(  A/100)  +  b  •  (A  -  100)  +  3.0  (2) 

where  a  and  b  are  coefficients  and  A  is  epicentral  distance. 

With  no  station  correction,  this  definition  gives  Mi  =  3.0  for  A  —  1  mm 
at  a  distance  of  100  km.  Richter’s  original  Ao  (attenuation)  values,  developed 
from  a  very  small  data  base  of  Wood-Anderson  seismometer  recordings  from 
southern  California,  are  surprisingly  correct,  as  demonstrated  recently  by  Bakun 
and  Joyner  (1984)  for  central  California  and  by  Hutton  and  Boore  (1987)  for 
southern  California. 

In  developing  a  new  Mi  scale  for  a  different  region,  it  is  first  of  all  necessary 
to  develop  a  new  attenuation  relationship  for  the  area,  in  terms  of  the  coefficients 
a  and  b  in  equation  (2).  A  problem  in  this  respect  is  that  regional  differences 
between  attenuation  at  100  km  makes  it  very  questionable  to  tie  the  two  scales 
together  at  this  distance.  In  using  a  shorter  distance,  such  as  10-20  km,  we  then 
on  the  other  side  run  into  a  problem  caused  by  the  fact  that  Richter’s  Ao  values 
are  more  poorly  defined  there. 

One  way  to  solve  this  problem  has  been  suggested  by  Hutton  and  Boore 
(1987),  who  developed  new  attenuation  relations  for  southern  California  from 
a  very  large  data  base.  Using  their  Ao  values,  which  are  identical  to  Richter’s 
at  a  distance  of  100  km,  a  magnitude  3.0  earthquake  should  be  recorded  with 
an  amplitude  of  10  mm  at  a  reference  distance  (Arf/)  of  17  km.  This  distance 
is  more  satisfactory  for  anchoring  purposes,  under  the  condition  that  our  local 
attenuation  function  can  be  evaluated  with  a  sufficient  precision  down  to  that 
distance.  If  that  is  not  the  case,  a  larger  A rej  must  be  chosen. 

This  leads  to  the  following  expression  for  Ao'. 

-  logAo  =  a-log(A/ Are/)  +  6-  (A  -  Are/)  +  K(AreJ)  (3) 

where  K(Are/)  is  determined  from  the  Ao  values  of  Hutton  and  Boore  (1987); 


28 


e.g.  A TCj  =  17  km  gives  K  =  2.0.  A  regression  analysis  should  then  be  aimed  at 
solving  for  the  parameters  a,  6,  5  and  Mi  in  equations  (1)  and  (3). 

The  collected  data  will  be  read  from  the  data  base  and  checked  manually  for 
noise  level,  clipping,  etc.  Fig.  2.3.3  shows  in  this  respect  a  panel  of  records  from 
an  Mi  3.1  earthquake  on  17  September,  1988,  with  epicentral  distances  from  176 
to  1308  km.  The  stations  include  three  from  SEISNOR,  NORESS  and  two  from 
VVNSN,  and  serves  as  a  good  illustration  of  the  usefulnes  in  mixing  records  from 
different  recording  systems. 

The  usable  records  will  then  be  corrected  for  instrument  response  (in  fre¬ 
quency  domain)  in  order  to  get  true  ground  motion  spectra  as  shown  in  Fig. 
2.3.4,  or  ground  motion  time  series,  in  the  cases  when  those  are  needed.  Simu¬ 
lated  Wood-Anderson  recordings  arc  then  established  by  applying  an  appropriate 
high  pass  filter  (2-pole  hp  at  1.25  Hz,  most  conveniently  applied  also  in  frequency 
domain),  and  the  resulting  time  series  is  then  plotted  (on  the  screen)  for  manual 
picking  of  maximum  amplitude.  These  amplitudes  are  then  used  regressively  by 
combining  equations  (1)  and  (3)  to  give 

N,  Ne 

logAtJ  =  -«-/ofif(A,y/Arf/)-<>-(A,y-Are/)-^5*:4j  +  X]A/i,^,(-A'(ArP/)  (4) 

k=  1  1=1 

where 

AtJ  =  simulated  Wood-Anderson  amplitude  of  earthquake  i  at  station  j 
=  Kronccker’s  symbol  (  =  1  if  i=j,  otherwise  0) 

N,  =  number  of  stations 

Ne  =  number  of  events 


The  parameters  to  be  determined  regressively  are  a,  6,  Sk  and  Mu  represent¬ 
ing  the  geometrical  spreading,  attenuation,  station  correction  and  magnitude, 
respectively. 

The  true  ground  motion  time  series  and  associated  Fourier  spectra  obtained 
from  the  data  base  serve  as  important  subsets  of  data  suited  for  further  wave 
attenuation  research.  Such  efforts  should  emphasize  the  interrelation  of  geomet¬ 
rical  spreading  and  anelastic  attenuation  in  the  computation  of  wave  attenuation, 
possible  local  and  regional  differences  in  wave  attenuation,  and  also  possible  az¬ 
imuthal  effects.  In  contrast  to  the  attenuation  terms  in  the  Mi  inversion,  the 
terms  here  can  be  evaluated  as  functions  of  frequency  (Dahle  et  al ,  1989). 
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Preliminary  results 

In  the  section  above,  we  have  outlined  a  procedure  for  simultaneous  inversion 
of  YVood-Anderson  amplitudes  vs.  epicentral  distance  in  order  to  estimate  A//, 
magnitudes,  station  corrections,  geometrical  spreading  and  anelastic  attenuation. 
This  can  be  done  at  the  same  time  or  through  a  multi-step  regression  analysis. 
In  order  to  test  the  stability  of  our  planned  procedure  we  selected  a  few  repre¬ 
sentative  recordings  as  shown  in  Table  2.3.1,  including  8  events  and  35  records. 
The  reference  distance  in  this  case  has  been  100  km,  the  estimated  magnitudes 
are  given  in  Table  2.3.1,  while  the  station  corrections  are  shown  in  Table  2.3.2. 

The  coefficients  a  and  b  in  equation  (4)  are  as  follow's: 

a  (geometrical  spreading)  :  0.70 

b  (anelastic  attenuation)  :  0.0010 

These  results  are  very  encouraging  in  that  we  for  all  of  the  21  estimated 
coefficients  get  reasonable  values  using  only  35  independent  observations.  In  gen¬ 
eral,  however,  we  must  expect  this  kind  of  analysis  to  be  more  unstable  than 
shown  hero,  and  that  therefore  the  various  restrictions  for  the  different  parame¬ 
ters  should  be  considered  closely,  in  addition  to  the  possibilities  for  a  multi-step 
regression  analysis.  The  resulting  magnitude  formula  would  read  as  follows  (ig¬ 
noring  the  station  corrections): 

Ml  =  log  A  +  0.70  •  log  A  +  0.001  •  A  +  1.5  (5) 

We  emphasize  that  these  results  are  preliminary  and  intended  only  for  test 
purposes. 

A.  Alsaker 
L.B.  Kvamine 
A.  Dahle 
H.  Bungum 
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Earthquake 

No.  of 
obs. 

Old 

Mi 

New 

Mi 

1986  Feb  5 

1 

4.8 

4.92 

1986  Oct  26 

1 

4.4 

4.44 

1987  Oct  31 

3 

3.6 

3.82 

1987  Nov  1 

6 

3.1 

3.18 

1987  Nov  3 

5 

2.3 

2.39 

1988  Jan  23 

10 

3.4 

3.35 

1988  Jan  31 

6 

3.8 

3.65 

1988  Aug  8 

3 

5.2 

5.20 

Table  2.3.1.  Earthquakes  used  in  testing  inversion  procedure  for  Mi  determi¬ 
nations,  with  number  of  records  for  each  event,  old  Mi  (from  current  formula), 
and  new  Mi  (from  inversion). 


Station 

No.  of 
obs. 

Correction 

NRS 

8 

-  0.02 

MORI 

1 

0.12 

FRS 

5 

-  0.02 

MOL 

4 

-  0.09 

LOF 

2 

0.17 

TRO 

4 

0.20 

KTK1 

3 

-  0.18 

SUE 

2 

0.05 

HYA 

2 

-  0.12 

ODDI 

2 

-  0.11 

KMY 

2 

-  0.02 

Table  2.3.2.  Station  corrections  from  testing  inversion  procedure  for  Mi  de¬ 
terminations,  with  station  name  in  column  1  and  number  of  records  in  column 
■» 
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Fig.  2.3.3.  Recordings  at  six  stations  from  an  Me  3.1  earthquake  on  17  Septcm 
ber  1988,  located  at  61.42°N,  1.57°E.  The  two  closest  stations  are  from  VVNSN, 
while  the  others  are  from  N0RESS  (NRS),  and  from  SEISNOR.  The  epicentral 
distances  are  SUE  17G,  HYA  249,  MOL  338,  NRS  543,  LOF  937  and  KTK  1308 
km,  respectively. 
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Fig.  2.3.4.  Fourier  acceleration  ground  motion  Lg  spectra  for  the  recordings  at 
the  stations  MOL  and  KTK  in  Fig.  2.3.3,  plotted  together  with  similar  estimates 
from  the  preceding  noise.  It  is  seen  that  KTK  (recording  an  Mi  3.1  earthquake 
at  a  distance  of  1308  km)  has  a  positive  signal-to-noise  ratio  only  around  2-G  Hz. 


2.4  A  generalized  beamforming  approach  to  real  time  network 
detection  and  phase  association 

Introduction 

The  objective  of  this  study  is  to  investigate  methods  for  multi-array  detection, 
phase  association  and  location,  using  as  primary  data  sources  the  three  regional 
arrays  (NORESS,  ARCESS,  FINESA)  in  Fennoscandia.  The  emphasis  is  on 
approaches  which  are  suitable  for  incorporation  in  an  expert  system  environment. 
As  part  of  this  effort,  the  problem  or  continuously  monitoring  the  regional  seismic 
noise  field  is  also  addressed,  with  the  purpose  to  obtain  a  quantitative  assessment 
of  the  upper  limit  of  magnitudes  of  seismic  events  that  would  go  undetected  by 
such  a  network. 

A  more  exhaustive  treatment  of  the  subject  of  this  contribution  is  offered  in 
our  Scientific  Report  No.  1  under  this  contract  (GL  Report  Number:  GL-TR-89- 
0171),  which  will  also  appear  in  the  December  1989  volume  of  the  Bulletin  of  the 
Srismnlogien!  Society  of  America. 

General  approach 

In  the  processing  of  seismic  network  data,  individual  phase  detections  corre¬ 
sponding  to  the  same  seismic  event  must  be  properly  associated  and  grouped  to¬ 
gether.  This  is  today  usually  done  starting  with  an  initial  trial  epicenter  and  then 
applying  various  search  strategies  supplemented  by  combinational  techniques. 

This  paper  presents,  and  gives  examples  of  application  of,  a  method  for  as¬ 
sociating  phase  detections  from  a  network  of  stations,  which  is  analogous  to  the 
conventional  delay-and-sum  beamforming  commonly  applied  in  array  detection 
processing.  A  number  of  beams  are  steered  to  a  predefined  grid  of  aiming  points 
in  a  geographical  space.  Each  beam  has  an  associated  set  of  time  delays,  where 
each  delay  corresponds  to  the  predicted  travel  time  for  a  given  phase  at  a  given 
station. 

We  assume  that  the  data  of  each  network  station  is  initially  subjected  to  a 
detection  processing  procedure,  whereby  a  list  of  phase  detections  and  attributes 
is  generated.  The  beamforming  process,  for  a  given  beam,  at  time  1',  can  be 
described  as  looking  for  a  pattern  of  detections/non-detections  that  matches  the 
predicted  pattern  for  a  hypothetical  event  with  origin  time  T  and  location  within 
the  beam  region.  The  actual  beam  value  is  derived  from  probabilistic  consid¬ 
erations,  and  in  essence  describes  how  well  the  observed  pattern  matches  the 
prediction.  By  moving  along  the  time  axis,  we  thus  obtain  a  beam  trace  that  can 
be  subjected  to  standard  threshold  algorithms  for  detection.  The  process  can  be 
supplemented  by  various  individual  “quality  of  fit”  measures  calculated  at  each 
time  point. 
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After  processing  individual  beams  as  described  above,  a  grouping/reduction 
process  is  applied  to  the  set  of  all  beam  detections  to  eliminate  side  lobe  detec¬ 
tions.  This  then  results  in  an  event  list,  comprising  origin  times,  locations  and 
a  list  of  associated  phases.  Further  refinement  can  then  be  achieved  by  standard 
techniques  for  accurate  hypocenter  determination,  magnitude  computation,  final 
consistency  check,  etc. 

The  generalized  beamforming  approach  also  provides  a  convenient  tool  to 
monitor  continuously  the  seismic  noise  field  associated  with  a  given  beam.  An 
application  of  particular  interest  in  a  monitoring  situation  would  be  to  calculate, 
at  each  step  in  time,  upper  confidence  limits  for  the  magnitude  of  non-detected 
events  for  each  beam.  This  would  be  useful  to  obtain  a  realistic  assessment  of 
actual  network  detection  capabilities,  at  any  given  point  in  time.  The  paper 
presents  an  example  of  practical  application  of  this  approach. 

Regional  phase  association 

The  method  has  been  applied  to  a  data  base  comprising  24  hours  of  record¬ 
ings  from  the  regional  arrays  NORESS,  ARCESS  and  FINESA,  with  a  beam 
deployment  covering  Fennoscandia  and  adjacent  areas. 

A  RONAPP-type  detector  was  first  applied  to  each  array  individually,  using 
the  broad  band  F-K  method  to  obtain  phase  velocity  and  azimuth  for  each  de¬ 
tected  phase.  The  resulting  detection  lists  then  provided  the  input  to  the  network 
processor. 

The  beam  grid  used  for  network  processing  is  shown  in  Figure  2.4.1,  and 
comprises  altogether  121  aiming  points,  approximately  equally  spaced.  Typical 
distance  between  aiming  points  is  150  km. 

In  the  network  beamforming  process,  a  simple  model  of  assigning  0/1  proba¬ 
bilities  to  individual  phases  at  each  station  was  used.  We  required  that  estimated 
phase  velocities,  azimuth,  dominant  frequency  and  arrival  times  fall  within  pre¬ 
defined  ranges  for  a  phase  detection  to  be  accepted  for  a  given  beam.  These 
tolerance  ranges  are  specified  in  Table  2.4.1.  Note  in  particular  that  only  very 
general  criteria  are  applied,  and  we  have  made  no  attempt  to  optimize  perfor¬ 
mance  by  regionalization. 

With  this  simplified  model,  the  network  beamforming  process  in  practice  was 
reduced  to,  for  each  beam  and  each  time  T,  counting  the  number  of  phase  matches 
for  a  hypothetical  event  located  in  the  beam  region  and  having  origin  time  T. 
The  detection  threshold  was  set  equal  to  2.  Thus,  all  occurrences  of  two  or  more 
matching  phase  detections  were  flagged  as  potential  events.  A  typical  beam  trace 
is  shown  in  Figure  2.4.2. 

A  grouping  procedure  was  then  applied  to  the  overall  beam  detection  list. 
This  was  done  by  successively  linking  together  entries  in  the  beam  detection  list 


in  such  a  way  that  a  new  entry  would  be  linked  if  it  had  at  least  one  individual 
phase  detection  in  common  with  a  previous  entry  in  the  group.  The  maximum 
allowable  duration  of  a  group  was  set  to  10  minutes  (in  practice,  the  longest 
duration  was  7  minutes  for  this  data  set).  In  order  to  resolve  obvious  multiple 
events,  groups  were  split  up  if  two  P-detections  from  the  same  array  occurred 
with  more  than  30  seconds  arrival  time  difference. 

The  results  are  summarized  in  Table  2.4.2.  It  is  important  to  note  that  the 
total  of  91  groups  comprise  all  possible  events  that  could  be  associated,  given 
the  station  detection  lists.  Also,  a  scrutiny  of  the  data  shows  that  only  3  of 
these  groups  contain  multiple  events,  all  of  these  being  small  presumed  mining 
explosions  semi  by  one  array  only. 

Some  of  the  entries  in  Table  2.4.2,  e.g.,  those  generated  from  two  secondary 
phases,  are  probably  questionable  seismic  events,  and  even  if  real,  may  be  im¬ 
possible  to  locate  accurately.  An  upper  magnitude  limit  could  be  estimated  for 
such  events,  in  order  to  determine  whether  further  detailed  analysis  is  desirable. 
However,  the  large  majority  of  the  entries  appear  to  correspond  to  real  seismic 
events,  and  the  grouping  procedure  facilitates  the  subsequent  detailed  analysis  of 
the  associated  phases. 

The  network  beamforming  procedure  gives  an  initial  cstimateof  event  location 
by  selecting  the  “best  beam"  in  each  group.  This  is  defined  as  the  beam  with 
the  greatest  number  of  associated  phase  detections,  and  if  equality,  the  smallest 
average  time  residual  of  the  detected  phases.  Since  the  initial  beam  grid  is  very 
coarse,  we  applied  a  beampacking  algorithm  for  each  detection  group,  using  a 
grid  spacing  of  20  km  in  order  to  improve  the  location  estimate. 

The  results  of  this  automatic  procedure  are  shown  in  Table  2.4.3  for  those 
events  for  which  independent  location  estimates  were  available.  We  note  that  the 
estimates  are  very  consistent  (median  difference  40  km),  and  thus  the  beam  results 
can  be  used  as  a  reasonable  first  estimate  of  event  location.  For  more  accurate 
results,  available  techniques  for  accurate  hypocenter  location  (e.g.,  TTAZLOC) 
should  be  used. 

Continuous  monitoring  of  upper  event  magnitude  limits 

As  a  second  application  of  the  generalized  beamforming  procedure,  we  now 
address  the  problem  of  monitoring  the  noise  levels  on  each  beam,  and  use  this 
information  to  assess  the  size  of  events  that  might  go  undetected. 

In  formulating  the  approach,  we  consider  a  given  geographical  location,  and 
a  given  “origin  time”  of  a  hypothetical  event.  Assume  that  N  seismic  phases  are 
considered  (there  might  be  several  stations  and  several  phases  per  station). 

For  each  phase,  we  have  an  estimate  5,  of  the  signal  (or  noise)  level  at  the 
predicted  arrival  time.  For  P-phases,  5,  might  be  the  maximum  STA  value  (1  sec- 
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one!  integration  window)  within  ±  5  seconds  of  the  predicted  time.  For  Lg,  .S', 
might  lie  the  average  STA  value  over  a  10-20  seconds  window. 

We  assume  that  the  network  has  been  calibrated  (or  alternatively  that  stan¬ 
dard  attenuation  values  are  available),  so  that  magnitude  correction  factors  (6,) 
are  available  for  all  phases.  Thus,  if  a  detectable  signal  is  present: 

m,  =  log(Si)  +  b,  (i  =  1, 2, . . .  N ) 

Here,  in,  are  estimates  of  the  event  magnitude  m.  Statistically,  we  can  consider 
each  m,  as  sampled  from  a  normal  distribution  (m,r).  (A  standard  value  of  r  = 
0.2  seems  reasonable  for  a  small  opicentral  area.) 

Let  us  now  assume  a  “noise  situation”,  i.e.,  that  there  are  no  phase  detections 
corresponding  to  events  at  the  given  location  for  the  given  origin  time. 

Wo  then  have  a  set  of  “noise”  observations  a,,  where 

a,  —  log( S  -  i)  +  b,  ( i  =  1, 2, . . .  A' ) 

If  a  hypothetical  event  of  magnitude  m  were  present,  it  would  have  phase 
magnitudes  m,  normally  distributed  around  m.  We  know  that  for  each  phase, 

TTlj  <  C:  (t  =  1,2, . . .  N) 

Let  us  look  at  the  function 

f(m)  —  Prob(all  m,  <  a./event  magnitude  m) 

For  each  phase 

f i ( m )  =  Prob( in,  <  a,/m)  =  1  -  4>( - — ) 

<T 

where  <1>  is  the  standard  (0,1)  normal  distribution. 

Thus,  assuming  independence, 

/(">) = n  /«(m) 

t= i 

The  90  per  cent  upper  limit  is  then  defined  as  the  solution  of  the  equation  f(m) 
=  1  -  0.90 

It  is  important  to  interpret  the  90  per  cent  limit  defined  above  in  the  proper 
way.  Thus,  it  should  not  be  considered  as  a  90  per  cent  network  detection  thresh¬ 
old  since  we  have  made  no  allowance  for  a  signal-to-noise  ratio  which  would  be 
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required  in  order  to  detect  an  event,  given  the  noise  levels.  Rather,  the  computed 
level  is  tied  to  the  actually  observed  noise  values,  and  to  the  fart  that  any  hypo¬ 
thetical  signal  must  lie  below  these  values.  Our  90  per  cent  limit  represents  the 
largest  magnitude  of  a  possible  hidden  event,  in  the  sense  that  above  this  limit, 
there  is  at  least  a  90  per  cent  probability  that  one  or  more  of  the  observed  noise 
values  would  be  exceeded  by  the  signals  of  such  an  event. 

As  an  application  of  the  method,  we  selected  an  area  as  shown  in  Figure  2.4.3 
situated  at  similar  distance  from  the  three  arrays.  For  each  of  the  three  arrays, 
one  Pn  beam  and  one  Lg  beam  were  steered  to  this  location.  The  beam  traces 
were  filtered  using  the  frequency  bands  :i  'i  Hz  (Pn)  and  2-1  Hz  (Lg).  Magnitude 
calibration  values  (/it)  were  obtained  by  processing  previously  recorded  events  of 
known  magnitude  (Ml)  and  at  similar  distance  ranges,  and  then  determining  5, 
values  independently  for  Pn  and  Lg. 

Based  on  these  input  traces  from  the  three  arrays,  a  network  beam  was  then 
formed,  using  time  delays  for  each  phase  that  corresponded  to  the  given  location. 
Arrival  time  tolerances  were  set  to  ±  5  seconds  for  P-phases  and  ±  10  seconds 
for  secondary  phases.  This  is  roughly  consistent  with  a  beam  radius  of  50  km  as 
shown  on  the  figure. 

We  chose  to  analyze  a  3  1/2  hour  interval  during  which  four  regional  seismic 
events  of  Ml  >  2.0  were  reported  in  the  Helsinki  bulletin.  No  events  were  re¬ 
ported  near  the  beam  area  in  this  period.  Figure  2.4.4  shows  90  per  cent  upper 
magnitude  limits  as  previously  defined,  plotted  as  a  function  of  time.  In  this  fig¬ 
ure,  only  the  Pn  phase  has  been  used,  and  the  three  arrays  are  shown  individually 
and  in  combination  (bottom  trace). 

It  is  clear  from  Figure  2.4.4  that  when  considering  individual  arrays  only, 
there  are  several  possible  time  intervals  when  relatively  large  events  {Ml  ~  2.0 
3.0)  might  go  undetected  because  of  signals  from  interfering  events.  However, 
when  the  Pn  phases  are  combined,  these  instances  occur  much  more  seldom. 

Figure  2.4.5  shows  a  similar  plot,  but  this  time  including  both  the  Pn  and  the 
Lg  phase  for  each  array.  Even  on  an  individual  array  basis,  this  causes  substantial 
reduction  in  the  upper  magnitude  limits.  For  the  combined  plot  (bottom  trace 
of  Figure  2.4.5),  which  takes  into  account  all  6  Pn  and  Lg  phases  from  the  three 
arrays,  we  see  that  the  upper  limit  is  well  below  Ml  =  2.0  for  the  entire  time 
interval.  Thus,  we  may  conclude  that,  at  the  specified  level  of  confidence,  no 
event  of  Ml  =  2.0  or  higher  occurred  in  the  beam  region  during  the  time  period 
considered. 

Conclusions 

With  regard  to  phase  association,  the  generalized  beamforming  technique 
provides  an  effective  method  to  group  all  combinations  of  individual  phase  de- 


41 


tections  that  could  possibly  correspond  to  the  same  seismic  event.  At  the  same 
time,  preliminary  estimates  of  epicenter  and  origin  time  are  obtained. 

'The  primary  importance  of  this  would  be  to  obtain  a  starting  point  for  sub 
sequent  detailed  interactive  analysis  aimed  at  precise  determination  of  source 
parameters.  In  particular,  expert  system  approaches  (either  script-based  or  rule- 
based)  could  be  invoked  at  this  stage.  The  advantage  of  applying  the  generalized 
beamforming  as  the  first  step  is  to  reduce  the  amount  of  combinational  process¬ 
ing  that  would  be  necessary  otherwise.  It  is  here  noteworthy  that  the  processing 
load  when  applying  generalized  beamforming  increases  in  a  linear  fashion  when 
the  number  of  individual  phase  detections  increase,  whereas  combinational  pos¬ 
sibilities  tend  to  increase  exponentially.  While  we  have  in  this  paper  used  only  a 
three-array  network,  the  extension  to  larger  networks  is  clearly  straightforward. 

The  application  of  the  method  to  provide  continuous  monitoring  of  upper 
magnitude  limits  at  specified  beam  locations  provides  a  useful  supplement  to 
standard  statistical  network  capability  studies.  In  particular,  this  application 
would  give  a  way  to  assess  the  possible  magnitude  of  non-detected  events  during 
the  coda  of  large  earthquakes.  In  such  situations,  it  would  be  appropriate  to  use 
global  network  data  and  include  as  many  relevant  phases  as  possible  for  each 
network  station.  For  example,  while  an  expected  P  phase  at  a  given  station  may 
be  obscured  by  the  earthquake  coda,  later  phases  such  as  PcP  or  PP  may  be  less 
influenced,  and  the  noise  level  at  their  respective  expected  arrival  times  would 
therefore  provide  important  information  as  to  the  size  of  possible  undetected 
events. 

As  a  final  comment,  we  note  that  the  approach  presented  here  to  upper  limit 
magnitude  calculation  could  be  applied  to  extend  the  utililty  of  various  discrim¬ 
inants,  such  as  At,  :  !/)(,.  For  small  explosions,  surface  waves  frequently  are  too 
weak  to  be  observed  at  any  station  of  the  recording  network.  Obtaining  reliable 
upper  bound  on  M,  in  such  cases  would  expand  the  Tange  of  usefulness  of  this 
discriminant.  In  practice,  an  “upper  bound”  for  single-station  measurements  has 
often  been  given  as  the  “noise  magnitude”  at  that  station,  i.e..  the  M,  value  that 
corresponds  to  the  actually  observed  noise  level  at  the  expected  time  of  Rayleigh 
wave  arrival.  The  proposed  procedure  will  include  this  as  a  special  case  of  a  more 
general  network  formulation. 

In  future  studies,  we  plan  to  investigate  the  application  of  more  sophisticated 
probabilistic  models  in  the  generation  of  beam  traces  and  the  continuous  ex¬ 
traction  of  features  associated  with  the  individual  beams.  Application  to  larger 
networks,  including  teleseismic  monitoring  using  global  network  data,  will  also 
be  considered. 


Frode  Ringdal 
Tormod  Kvaerna 
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Pn 

Phase  Type 
Pg  Sn 

Lg 

R-g 

Distance  interval1) 

(km)  for  which  a  phase 
is  accepted 

160-3000 

0-600 

160-3000 

0-2000 

0-400 

Maximum  allowable 
deviation  from  predicted 
arrival  time  (s) 

15 

20 

30 

35 

40 

Maximum  allowable 
azimuth  deviation 
(degrees) 

20 

20 

20 

20 

20 

Acceptance  limits  for 
apparent  phase 
velocity  (km/s) 

5.8- 1-1 

5.8-10 

3. 2-5.8 

3. 0-5.0 

2. 5-3. 7 

Acceptance  limits2) 
dominant  frequency 

(Hz) 

0.5-20 

0.5-20 

0.5-20 

0.5-20 

0.5-20 

*)  For  NOR  ESS,  the  Rg  phase  is  not  included  in  the  phase  table 
2)  For  FINESA,  a  lower  frequency  limit  of  0.9  Hz  is  used  for  all  phases. 


Table  2.4.1.  Acceptance  limits  for  parameters  used  in  the  network  beamforrning 
process. 
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— _ _ 1 

Number  of  phases  for 
best  beam  in  each  group 

INo.  of  phase  groups:  | 

2 

3 

4 

5 

6 

7 

8 

NORESS  only 

18 

13 

4 

1 

0 

0 

0 

0 

ARCESS  only 

34 

19 

10 

4 

1 

0 

0 

0 

E 'IN ESA  only 

14 

13 

1 

0 

0 

0 

0 

0 

Two  arrays 

17 

9 

4 

3 

0 

1 

0 

0 

Three  arrays 

8 

0 

0 

2 

0 

1 

3 

2 

Totals 

91 

54 

19 

10 

1 

2 

3 

2 

Table  2.4.2  Phase  groups  associated  by  the  network  beamforming  procedures 
for  a  24-hour  interval. 
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Event 

No. 

Date 

Time 

Network 
Lat.  Lon. 

Mag. 

Ml 

No.  of 
phases 

No.  of 

arrays 

Beamforming 
Lat.  Lon. 

‘Error* 

(kin) 

1 

88/03/17 

08.40.25.0 

57.73 

11.03 

2.5 

7 

3 

57.9 

10  4 

39 

2 

08.46.18.7 

58.07 

11.36 

2.6 

6 

2 

57.9 

10.8 

:5(. 

3 

09.07.10.3 

58.08 

11.43 

2.7 

8 

3 

57.8 

10.8 

47 

4 

10.21.23.0 

69.60 

29.90 

2.9 

8 

3 

69.6 

30.5 

5 

10.27.20.0 

59.20 

27.60 

2.3 

4 

2 

59.5 

27.5 

34 

0 

10.46.21.0 

59.20 

27.60 

<2.0 

2 

1 

59.8 

28.7 

99 

7 

11.18.48.0 

59  30 

27.20 

2.3 

5 

3 

58.9 

26.7 

•  j’J 

8 

11.54.41.0 

65.80 

24.70 

<2.0 

5 

1 

66.6 

24.4 

S‘J 

9 

11.57.57.9 

60.57 

8.36 

1.8 

2 

1 

60.6 

8.1 

14 

10 

12.02.30.0 

59.40 

28.50 

2.1 

3 

2 

59.5 

2S.2 

'JO 

11 

12.42.22  9 

59.78 

10.76 

2  3 

3 

1 

59.5 

10.0 

•VJ 

12 

14.13.14.0 

58.33 
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5.57 

3.2 
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3 

60.0 
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36 

Table  2.4.3.  Location  estimates  obtained  automatically  from  the  beam  packing 
procedure  compared  with  independent  network  locations  from  the  Helsinki  and 
Bergen  bulletins.  Note  the  good  consistency,  especially  for  events  with  more  than 
one  detecting  array. 
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LATITUDE  (DEG  N) 


Fig.  2.4.1.  Roam  grid  used  in  the  generalized  beamforining  procedure  for  the 
purpose  of  associating  regional  phases  from  NORESS,  ARCESS  and  F1NESA. 
The  location  of  the  three  arrays  is  shown  on  the  map. 
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BEAM  LOCATION:  57.0N  10.0E 


Fig.  2.4.2.  Example  of  typical  output  trace  for  one  network  beam  (steered  to 
57"N,  10°E).  In  the  3  1/2  hour  interval  shown,  there  were  4  confirmed  seismic 
events  located  in  the  beam  region.  These  were  all  correctly  detected  (arrows), 
and  no  false  phase  associations  occurred  during  the  interval. 
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LATITUDE  (DEG  N) 


CONTINUOUS  THRESHOLD  MONITORING 


Fig.  2. 4. 3.  Location  of  the  '-cam  area  used  in  the  example  of  continuous  moni¬ 
toring  of  upper  magnitude  limits  on  non-detected  events.  The  area  covers  a  circle 
of  approximately  50  km  radius,  and  is  situated  at  similar  distances  om  the  three 
arrays. 
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CONTINUOUS  THRESHOLD  MONITORING  -  PN  PHASE 


90%  PROBABILITY 


oaoo  0900  10.00  it  oo 

hour 


Fig.  2.4.4.  Results  from  the  continuous  threshold  monitoring  of  the  area  shown 
in  Fig.  2.4.3  for  a  3  1/2  hour  period,  using  Pn  phases  only.  The  top  three  traces 
show,  for  each  array,  the  largest  magnitude  of  a  possible  non-detected  event 
(confidence  90  per  cent)  as  a  function  of  time.  The  bottom  trace  shows  the  result 
of  combining  the  observations  from  all  three  arrays  (Pn  phase  only)  as  described 
in  the  text. 


CONTINUOUS  THRESHOLD  MONITORING  -  PN  AND  LG  PHASES 


08.00  09.00  10.00  11.00 
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Fig.  2.4.5.  Same  as  tig.  ‘2.4.4,  but  using  both  the  Pn  and  Lg  phases  for  the 
upper  magnitude  limit  calculations.  Comparing  with  Fig.  2.4.4,  we  note  that  this 
serves  to  lower  the  thresholds,  both  for  each  individual  array  (top  three  traces) 
and  for  the  combined  results  (bottom  trace). 
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2.5  Region-specific  knowledge  derived  from  analysis  of  NORLSS 
data 

Introduction 

rite  Intelligent  Monitoring  System  (IMS)  that  is  currently  being  installed  at 
N'ORSAR,  provides  functions  that  allow  specific  region-dependent  knowledge  to 
lie  taken  into  account  in  the  automatic  processing  of  multi  array  data.  To  provide 
such  knowledge  for  enhancement  of  the  performance  of  IMS  i>  and  will  remain 
an  important  aspect  of  the  reseat  ch  at  N'ORSAR  under  this  <  out  ran.  In  Miis 
colli  I  i  1)11 1  ion  ,  we  lies!  l  ilie  t  wo  pro  jects  that  are  aimed  at  deliving  leg  II  III  s  pe<  i  lie 
knowledge  for  the  benefit  ol  IMS. 

Statistics  on  nt/hinal  phase  attribuh  s  from  \UIIESS  detection  lists 

The  NORKSS  detection  lists  for  the  period  19S5  to  1988  contain  more  than 
200,000  entries,  comprising  leleseismic,  regional  and  local  phase  arrivals.  Kite 
network  bulletins  for  this  four-year  period  published  by  the  University  of  Bergen, 
Norway,  and  the  University  of  Helsinki.  Finland,  contain  altogether  10.000  events 
that  are  local  or  regional  to  the  NO  RUSS  array.  For  each  of  these  1 G. 000  events, 
we  calculated  tlm  expected  arrival  times  at  NORKSS  for  the  regional  phases 
Fit.  Pg.  S'  and  l.g  and  then  searched  the  detection  lists  to  see  whether  the 
theoretical  arrival  times  matched  those  ol  actual  phase  arrivals  detected  on  the 
NORKSS  array.  Certain  criteria  were  established  in  this  regard,  and.  e.g..  lor 
the  Fn  phase  these  were  that  the  arrival  times  should  not  deviate  by  more  than 
(>  seconds  from  the  predicted  ones,  the  phase  velocity  (routinely  determined  by 
frequency- wavenumber  (f-k)  analysis  for  each  detected  signal)  should  be  within 
the  interval  0  15  ktn/s,  and  the  arrival  azimuths  (also  determined  by  f  k  analysis) 
should  be  within  10"  of  the  predicted  values. 

This  process  of  merging  of  the  network  bulletins  with  the  NORF.SS  detection 
lists  prodme.s  a  wealth  of  informal  ion  on  attributes  like  phase  velocities,  arrival 
time  residuals  and  arrival  azimuth  residuals  relative  to  the  network  locations. 
The  re.Milt.s  may  also  be  used  as  a  basis  for  compiling  data  bases  of  interesting 
even's  from  specific  source  regions. 

As  a  result  of  this  merging  process,  approximately  9.000  arrivals  detected  on 
the  NORKSS  array  could  be  classified  as  either  Pn,  Pg,  Sn  or  Lg  originating  from 
these  10,000  events.  Figure  2.5.1  (top)  shows  the  number  of  Pn  phases  detected 
at  NORKSS  from  these  events  and  how  the  corresponding  events  are  distributed 
on  a  grid  of  Tx  2"  blocks  (north-south  and  east-west,  respectively).  The  region 
around  Bergen  on  the  western  coast  of  Norway,  the  Kstonia -Leningrad  region  of 
the  western  USSR,  the  Finland- USSR  border  region  at  around  •t()‘T..  05" N.  an 
area  in  northern  Sweden,  and  the  Kola  Peninsula  of  the  USSR  stand  out  clearly 
in  this  map.  The  pattern  observed  partly  reflects  the  station  distribution  ol  the 


reporting  agencies,  hut  more  basically  gives  an  overview  of  the  main  mining  areas 
in  the  Ne, relic  countries  and  the  northwestern  USSR.  Figure  2.-r i.l  (bottom)  also 
shows  the  average  phase  velocity  for  Pn  phases  from  events  within  each  block.  We 
see  that  these  velocities  fall  within  the  expected  range  around  8  km/s  for  most 
regions  of  this  map.  An  exception  is  the  Estonia- Leningrad  area  of  the  western 
USSR,  where  Pn  phase  velocities  averaging  more  than  10  km/s  are  observed. 
Arrival  azimuth  residuals  for  the  Pn  phases  arc  shown  in  Figure  2.5.2.  It  is  seen 
that  the  high  Pn  phase  velocities  observed  for  the  Estonia- Lett  ingrad  region  are 
accompanied  by  a  rather  complex  pattern  of  azimuth  residuals. 

The  azimuth  residuals  for  the  I.g  phase  are  shown  in  Figure  2.5.3.  It  is  note- 
worths  that  the  Estonia-  Leningrad  region  exhibits  rather  moderate  I.g  azimuth 
residuals. 

It  should  be  noted  that  phase  velocities  and  arrival  azimuths  for  this  study 
were  derived  using  tint  conventional  or  narrow-band  f-k  analysis  method.  Since 
1989,  the  broad-band  f-k  method  is  used  in  the  analysis  of  all  array  data  at 
NORSAR.  and  will  also  be  used  by  the  IMS.  Several  studies  have  testified  that 
the  broad  band  approach  is  more  stable  than  the  narrow-band  method,  so  a  new 
investigation  along  the  lines  of  this  study  should  be  undertaken,  based  on  the 
detection  lists  for  1989. 

.4  detailed  study  of  lO.i  events  in  the  western  Norway  ana 

Another  investigation  aimed  at  deriving  region-specific  knowledge  from  anal¬ 
ysis  of  NOR  ESS  data  was  a  detailed  study  of  103  events  in  the  western  Norway 
a  ret  i .  The  locations  of  these  events  are  shown  in  Figure  2.5.4.  The  large  major¬ 
ity  of  the  events  are  in  the  distance  interval  250-700  km  and  azimuth  interval 
2 1 01'  330''.  relative  to  NORESS.  All  phases  that  could  be  picked  for  each  of  the 
103  events  were  subjected  to  f-k  analysis,  using  both  the  narrow  -band  and  broad¬ 
band  estimation  method.  The  arrival  azimuths  derived  from  these  analyses  were 
compared  with  at  rival  azimuths  calculated  from  the  network  locations  published 
by  the  University  of  Bergen.  Norway,  and  the  British  Geological  Survey. 

Arrival  azimuth  residuals  (defined  as  azimuth  estimated  from  f-k  analysis 
minus  azimuth  according  to  the  network  bulletin)  for  these  events  lor  the  Pit. 
Pg,  Sn  and  I.g  phases  are  plotted  as  function  of  network  azimuths  in  Figure;  2.5.5 
for  the  broad  band  estimation  method,  and  in  Figure*  2.5.0  for  the  narrow-band 
estimation  method.  It  is  dearly  seen  that  the  broad-band  method  provides  the 
more*  stable  estimates,  and  the  standard  deviation  of  the  azimuth  residuals  for 
the  broad- baud  method  were  found  to  be  of  the  order  of  5-7  degrees,  as  opposed 
to  7-11  ele’grees  for  the  narrow-band  method.  The  smallest  standard  deviations 
were-  found  for  tie*  crustal  phases  Pg  and  Lg. 

An  interesting  observation  is  that  for  the  Pn  and  Sn  phases,  tlu*  azimuths 
(estimated  by  the  broad-band  method)  deviate*  systematically  from  the*  ‘true* 


values,  in  the  sense  that  the  residuals  are  negative  for  azimuths  in  the  range  210"- 
270°  and  positive  for  azimuths  between  270°  and  .120".  These  results  indicate  t  hat 
laterally  varying  structures  near  the  crust-mantle  boundary  rather  than  structure 
within  the  crust  are  important  in  the  sense  that  they  influence  the  azimuths 
estimated  at  NORESS.  More  work  is  needed  in  order  to  gain  further  insight  and 
understanding  of  these  effects. 

Conclusions 

The  merging  of  the  complete  NORESS  detection  lists  for  the  period  I  ?).8.r>-  1  OSS 
with  the  regional  network  bulletins  results  in  a  large  amount  of  statist  ical  infor¬ 
mation  on  characteristics  of  regional  phases  in  Eennoscandia  and  adjacent  areas. 
Similar  information  will  result  from  merging  with  the  AIKT'.SS  and  I  INKS  A  de¬ 
tection  lists.  The  knowledge  gained  from  this  exercise  will  be  of  relevance  and 
importance  to  the  IMS,  and  it  is  anticipated  that  this  knowledge  will  be  repre¬ 
sented  in  IMS's  knowledge  base. 

A  special  study  of  lO-'f  events  recorded  at  NORESS  from  the  western  Norway 
area  is  another  example  of  what  kind  of  studies  that  must  be  undertaken  in  order 
to  obtain  region-specific  knowledge  needed  by  the  IMS.  Among  other  results,  this 
special  study  established  that  the  arrival  azimuth  residuals  observed  at  NORESS 
for  the  Pn  and  Sn  phases  vary  rather  systematically  with  source  region  in  the 
western  Norway  area. 
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Fig.  2.5.1.  Top:  Number  of  i’n  pli  uses  ( () |!S  )  deter! i>, 
various  source  areas,  liottom:  ( '01  irip.mding  average  Pi 
location  ol  \OKPSS  is  indicate. 1  l»v  a  special  symbol. 
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Fig.  2.5.4  .  I'lio  map  shows  the  location  of  the  103  events  in  the  western  Norway 
legion  used  in  this  study.  The  location  of  NORESS  is  indicated  by  a  special 
symbol. 
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Fig.  2.5.5.  1  lie  figure  shows  the  azimuth  residuals  as  function  of ‘true’  azimuths, 
for  the  broad-band  estimation  method,  i  lie  residuals  am  given  (from  top  to 
bottom)  for  the  l’n.  l’g,  Sn  and  Lg  phases. 
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Fig.  2.5.6.  Same  as  Fig.  2.5.5,  but  for  the  narrow-  band  estimation  method. 
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